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National Bureau of Standards 
Materials Measurements 

Summary 

Ibis report describes BBS work for NASA in support of NASA's Materials 
Processing in Space Program under NASA Government Order H-27954B (Properties 
of Electronic Materials) covering the period April 1, 1982 to March 31, 1983. 
This work is directed toward measurement of tutorials properties important 
to the design and interpretation of space processing experiments and 
determinations of how the space environment say offer a unique opportunity 
for performing improved measurements and producing materials with improved 
properties. 

The work is being carried out in three independent tasks: 

Task 1. Surface Tensions ana Their Variations with Temperature and 
Impurities 

Task 2. Convection During Unidirectional Solidification 

Task 3. Measurement of the High Temperature Tbermophyeical Properties 
of Tungsten Group Liquids and Solids 

The results obtained for each task are given in detailed summaries in the 
body of the report. Baphaeis in Tasks 1 and 2 la on how the reduced 
gravity obtained in apace flight can affect convection and solidification 
procasaea. Emphasis on Taak 3 is toward development of techniques for 
thermodynamic meaaurcaenta on reactive materials, requiring levitation 
and containerless processing. 

With the advent of the Space Shuttle, it may become feasible to 
exploit the unique microgravity environment of apace flight to produce 
improved materials and improved measurements of Important materials 
properties. In materials processing on earth, gravity frequently produces 
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iensity-dr iven convection, thereby causing liquids to be stirred as they 
solidify. This stirring disturbs the quiescent boundary layer at the 
solidifying interface, and can be very undesirable if nearly perfect 
crystals are required. For example, it creates interface instabilities, 
introduces segregation of components and produces crystal defects in the 
resulting solid material. These defects and inhomogeneities, which are 
particularly troublesome in electronic technology and other advanced 
technical applications, might be avoided in materials produced under 
microgravity conditions. 

In Task 1, surface tension and surface segregation effects are being 
investigated to determine the possible influence of these surface effects 
on convection and solidification processes. In particular, measurements 
are being made of the surface tension of liquid silicon. Silicon is one 
of the most important materials used in electronic technology. Because 
of the need for fully defect-free material, not readily obtainable under 
earth-bound conditions where gravity-driven convection can be important, 
silicon is a major candidate material for space flight processing. Under 
microgravity conditions, Marangoni convection arising fro© surface tension 
gradients could be the main source of fluid flows affecting solidification 
processes. Thus surface tension information will be Important in planning 
materials processing in space. Because of the high reactivity of silicon, 
i ts surface tension is poorly known. Moreover, despite the technological 
importance of silicon, the dependence of its surface tension on temperature 
and impurities was really not known at all prior to the current work. 

Ibis work is designed to provide these surface tension measureiwuts. 
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In Task 2, solutal convection* which results from simultaneous 


temperature and composition gradients, and other convective phenomena 
arising during directional solidification are being investigated both by 
theoretical calculations and by experiments. The experiments are designed 
to test the theories and provide quantitative data on convection processes. 
Direct measurements of convection are being made during unidirectional 
solidification of transparent materials, particularly succinonitrile 
containing small amounts of ethanol. The theoretical calculations 
specifically include determinations of the effects a reduction in 
gravitational force will have on convection. 

Interactions between temperature gradients and composition gradients 
can make it impossible to avoid solutal convection and other density- 
driven convection effects during materials processing on earth. Since 
this type of convection depetfds oq gravity, its effects will be much less 
pronounced in space. KBS work Is directed toward provision of measurements 
and development of models to define the conditions under which this type 
of convection will be Important and determine how it can be avoided or 
controlled. 

In Task 3, assistance is being provided to a joint project involving 
Investigators from Rice University (Prof. J. Margrave) and General Electric 
Co. (Dr. R. T. Frost) in which a General Electric electromagnetic levitation 
facility is being applied to develop levitatlon/drop calorimetry techniques 
and determine their possibilities for use in space flight experiments. 

The critical advantage of such levitation experiments is the avoidance of 
specimen contamination since in highly reactive materials even slight 
contact with a container during heating and melting could completely 
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invalidate the experimental results. Measurements of the specific heat ; 

of liquid tungsten uncontaminated by reaction with container walls 

i 

currently are being pursued in earth's gravity experiments to look at the • 

limitations imposed by gravity in such work. For example, in these 

experiments, the large electromagnetic forces required to counteract ^ 

gravity tend to produce instabilities in the liquid droplets. A number 
of the automated techniques being developed in this work should also be 
useful for space flight experiments. 
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Task 1 


Surface Tensions and Their Variations 
with Temperature and Impurities 

S. C. Hardy 
Metallurgy Division 
Center for Materials Science 

Introduction 

Surface tension gradients produced by temperature or concentration 
variations on free liquid surfaces result in shear stresses which can in 
many situations set the fluid into motion. In low gravity materials 
processing the primary flu'** flows are expected to be of this type because 
buoyancy is greatly reduced. These thermocapillary flows can influence 
the temperature and solute fields in the bulk liquid and the morphology, 
microstructure and homogeneity of a solid produced from this liquid. 

Thus, knowledge of surface tension gradients Is essential to the under- 
standing of materials produced from the processes involving free liquid 
surfaces in l<»w gravity. We are measuring the surface tension of liquid 
silicon in support of space experiments with this material vhich are 
being developed. The results of these measurements will also be of 
Interest for silicon crystal growth research and production on earth. 

The temperature and chemical concentration dependences of the surface 
tensions of many materials are poorly known; for highly reactive, high 
melting point elements like silicon often only estimates of these parameters 
can be found In the literature. Thus in a recent computer study of the 
fluid flow in Czochralski growth of silicon, calculations were performed 
for values from -0.1 mJ/m^K to -0.4 mJ/m^K for K, the rate of change of surface 
tension with temperature (1). The flew was found to vary significantly 
both in magnitude and distribution for this variation in K. These results 
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show that accurate surface tension data is ecsential for meaningful calcu- 
lations or quantitative interpretation of observed flow phenomena in 
floating zones of silicon. 

The uncertainty in K for silicon stems basically from the lack of 
systematic measurements of the surface tension over a temperature range. 
Various workers have measured the surface tension at different temperatures 
[2, 3, A]. However, K values calculated from their data are highly 
inaccurate because systematic errors in the different experiments are 
sufficiently large to obscure the surface tension differences due to 
temperature. It is essential in determinations of K to measure the 
surface tension over a temperature range in the same experiment so that 
the effects of systematic errors are minimized. This is especially true 
for silicon because of its extreme chemical reactivity. Liquid silicon 
dissolves almost all materials to some extent thereby picking up impuric 
which can surface segregate and significantly lower the surface tension. 
Thus it must be assumed that most experiments will have unique impurity 
concentrations. 

The only systematic study of the surface tension of silicon over a 
temperature range that we know of is that of Lukin, et al [5]. The 
value they found for K is -0.1 mJ/a^K. Although this is a reasonable 
value, no details or data are presented. Moreover the density variation 
with temperature used in their work is significantly different from that 
generally accepted. Thus it is not possible to make any Judgment as to 
the accuracy of their measurements. 

Experimental 

We are using the sessile drop technique to measure the surface 
tension of liquid silicon. This method is based on a comparison of the 
profile of a liquid drop with the profile calculated by solving the 
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Young-Laplace equation. The comparison can be made in several ways; we 
are using the traditional Bashf orth-Adams [6] procedure with modern drop 
shape tables we have calculated which virtually eliminate interpolation 
errors. 

Although vre have used the ee Mile drop technique previously to 
measure the surface tension of gallium [7], the present experiments 
require essentially new apparatus because of the high temperatures 
involved. The melting point of silicon is 1410°C and it is desirable to 
make measurements to at least 160O°C in order to determine the temperature 
dependence of the surface tension. 

Figure 1 is a schematic diagram of the experimental apparatus. The 
silicon is contained in a cup which stabilizes the liquid drop and imposes 
a circular cross sectior, on the drop base, thus reducing errors rising 
from anisotropic wetting often encountered using flat drop supports. The 
outside diameter of the cup is used as the length standard in the 
calculation of the surface tension. The cup and drop are at the center 
o* a cylindrical susceptor which is heated inductively. The quartz jacket 
enclosing the susceptor and drop connects through an O-ring coupling to a 
vacuum system consisting of a water cooled baffle and diffusion pump, 
gauges, valves for admitting gases to the system, and a forepump. Quartz 
windows at both ends of the tube allow simultaneous photography and 
temperature measurement of the drop with an automatic pyrometer. 

Figure 2 shows the central region of the heating tube in more detail. 
The cup containing the silicon sits on a tungsten platform which in turn 
rasts on two tungsten rods. These rode extend well out of the susceptor 
into a much cooler region o f ^ne heating tube. They are clamped to a 
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Young-Laplace equation* the compar. ^on can be made in several ways; we 
are using the traditional Bashforth-Adams [6] procedure with modern drop 
shape tables we have calculated which virtually eliminate interpolation 
errors. 

Although we have used the sessile drop technique previously to 
measure the surface tension of gallium [7], the present experiments 
require essentially new apparatus because of the high temperatures 
involved. The melting point of silicon is 1410°C and it is desirable to 
make measurements to at least 16O0°C in order to determine the temperature 
dependence of the surface tension. 

Figure 1 is a schematic diagram of the experimental apparatus. The 
silicon is contained in a cup which stabilises the liquid drop and imposes 
a circular cross section on the drop base, thus reducing errors arising 
from anisotropic wetting often encountered using flat drop supports. The 
outside diam e ter of the cup is used as the length standard in the 
calculation of the surface tension. The cup and drop are at the center 
of a cylindrical susceptor which is heated inductively. The quarts jacket 
enclosing the susceptor and drop connects through an O-ring coupling to a 
vacuum system consisting of a water cooled baffle and diffusion pump, 
gauges, valves for admitting gases to the system, and a forepuotp. Quarts 
windows at both ends of the tube allow simultaneous photography and 
temperature measurement of the drop with an automatic pyrometer. 

Figure 2 shows the central region of the heating tube in more detail. 
The cup containing the silicon sits on a tungsten platform which in turn 
rests on two tungsten rods. These rods extend well out of the susceptor 
into a much cooler region of the heating tube. They are clamped to a 
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stainless steel block at one end and rest freely on an Identical block at 
the other end. These steel blocks rest on the vails of the quartz tube. 

The platform and cup can be moved along the rods and the entire assembly 
translates inside the quartz tube. This device permits convenient and 
rapid sample insertion and removal. 

The induction heater operates at 450 kHz and has a pover rating 
of 2.5 kW. This has proved adequate for these experiments if the 
susceptor is insulated sufficiently. The most promising insulation we 
have tried is zirconla, Zr0 2 » because of its low thermal conductivity and 
outgassiog rate. However, the stability of the insulation in conjunction 
with the susceptor material is critical. In our initial measurements ve 
used zirconla wrapped about a tantalum susceptor. At operating temperature 
the tantalum reduced t v e zirconla and generated volatile species according 
to the following reaction [8]: 

Ta(s) + Zr0 2 (s) - TaO(g) + ZrO(g) 

Similar reactions can occur with other combinations of Insulator and 
susceptors (9). At present ve are using zirconia insulation with molybdenum 
susceptors. This combination is reported to be non-reacting up to 2200°C 
[9, 10]. Metallic heat shields cannot be used because the skin depth at 
450 kIL is only 0.06 cm. Since the field decays exponentially aa the 
thickness over the penetration depth, a heat shield 0.006 cm thick reduces 
the field by 10Z, a power loss which would be intolerable with our equipment. 

The suaceptor is constructed of two concentric cylinders of molybdenum 
with the third, inner cylinder of tantalum because of its lower vapor 
pressure. The total thickness is about 0.11 cm, which is sufficient to 
reduce the field to about 14Z of the maximum val le. Thus the interior of 
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the susceptor is not field free. This a ay result in net forces on the 
sessile drop which will cause distortions. Although we do not have 
evidence for such distortions, it is desirable to decrease this field 
penetration. 

Because of the very high reactivity of liquid silicon, the choice of 
cup material is a critical feature of these measurements. Extensive 
studies of the compatibility of liquid silicon with various materials 
have been made recently in conjunction with the solar cell program of the 
DOE. This research suggests that boron nitride is the most promising of 
a very limited group of materials [11]. Liquid silicon has a contact 
angle with boron nitride of about 130° [11]. For our purposes a high 
contact angle is desirable because it minimizes the wetting of the top 
surface of the cup by the silicon with the attendent distortion of shape. 
More importantly, silicon does not dissolve boron nitride. Thus the 
silicon is not grossly contaminated by cup material. The bonding between 
the silicon and the boron nitride is sufficiently weak to preserve the 
cup on cooling. Most materials bond strongly to liquid silicon and are 
cracked by the differential thermal contraction on cooling. Thus the 
boron nitride cups can be used for many experiments. Although the DOE 
studies shoved that silicon did acquire some free boron when in contact 
with boron nitride, we do uot believe this will have a major effect on 
our measurements because boron has a higher surface tension than silicon 
and consequently will not surface segregate. The development of a bulk 
boron concentration will only Increase the surface tension. Studies of 
the boron silicon system indicate that 1 atomic percent boron will raise 
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the surface tension of silicon about 25 mJ/m^, i.e., approximately 3Z (5). 
This rould be a very high bulk Impurity concentration. 

The silicon used in these early experiments is polycrystalline 
material of six nines purity. The samples are cut from a large boule 
using a carbide cut-off wheel. After rinsing in water and alcohol, the 
silicon paralleplped is melted and solidified on a flat boron nitride 
block to form a button with roughly the shape of a sessile drop. This 
button is then placed in a boron nitride cup and a sessile drop is formed. 
The temperature of the silicon is raised to the TMaximxm and the drop is 
photographed as it is cooled to various temperatures by reducing the 
power setting of the induction heater. 

Results and Discussion 

Figure 3 shows the surface tension of liquid silicon as a function 
of temperature for three different samples. Two sets of data are In good 
agreement with each other: the third set lies about 40 mJ/m^ higher than 

these points. The slopes, however, are nearly the same for the two 
groupings of data. Linear regression analysis gives slopes of 

-0.245 ± 0.012 and -0.237 t 0.025 mJ 
m^K m^K 

for the lower and higher data respectively where the range is the estimated 
probable error. Thus our early measurements indicate the variation with 
temperature of the surface tension of silicon is approximately -0.24 mJ/m^K. 
This data was acquired before the construction of the vacuum system; the 
measurements were made in flowing argon which had been purified by passage 
over hot titanium. The 40 mJ/m^ difference in the two sets of data is 
probably due to a different partial pressure of oxygen in the flowing argon. 
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It is possible, of course, that the surface tension difference arises 
rron a drop asymmetry which we cannot observe or to bulk contamination 
with boron. These values of surface tension and the temperature variation 
of the surface tension should be regarded as preliminary and as a 
developmental step rather than a quantitative measurement. 

We also show in Figure 3 as a dashed line the results of the Russian 
measurements referred to previously which have only been reported in 
summary form [5]. Ibis curve is in fair agreement in absolute value of 
the surface tenaio.i with our lower data. However, the temperature 
variation is -0.1044 aJ/m^K, a value less than half that which we found. 

Part of the discrepancy is due to the use of different density data. 

Hoe authors of reference [5] claim to have measured the density and 
report their results by giving a value for the coefficient Kp which 
occurs in the expression 

p ■ "mo * *» (T - :4l0) - 

The value of for pure silicon they find is 0.1291 x 10' ^gm/cm^K. 

We have used in our calculations of surface tension the density 
measurements of Lucas [12] et al which are reported in detail, are 
systematic, and in good agreement with less extensive measurements by 
others, Lucas finds K p - 0.35 x 10~^gm/cm^K, a value significantly 
diftviient from the Russian work. If the lower Kp Is used to calculate 
surface tensions with our drop measurements, the temperature variation of 
the surface tension la reduced to -0.185 mJ/n^, a result still significantly 
higher than that reported in the Russian work. 

In summary, we have developed the instrumentation required to measure 
the temperature variation of the surface tension of silicon. Preliminary 
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Measurements in flowing argon using boron nitride cups to contain the 
drops suggest that this temperature variation will be near -0.24 mJ/m^K. 

In future work we will be asking measurements in vacuum and in atmospheres 
containing controlled concentrations of hydrogen and oxygen in order to 
determine the sensitivity of the silicon surface tension to oxygen partial 
pressure. In addition, the effects of other container materials will be 
investigated. Of particular interest is Si02 because of its wide use in 
industrial silicon crystal growth. Some determination of the impurity 
content of the silicon after the surface tension measurements is highly 
desirable to assess the extent to which the silicon is contaminated by 
contact with the cup. Ws will characterise our samples using Auger 


spectroscopy. 
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Schematic diagram of experimental apparatus: (1) quartz tube, 

(2) induction coil, (3) O-ring aeal, (4) window. 
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Fig. 2. Schematic diagram of high temperature zone of experimental 

apparatus: (1) quartz tube, (2) zirconia insulation, (3) molybdenum 

susceptor, (4) W rods, (5) W tray, (6) BN cup. 
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The surface tension of silicon a a a function of temperature 
(preliminary results). 






Task 2 


Convection During Unidirectional Solidification 

R. J. Schaefer and S. R. Coriell 
Metallurgy Division 
Center for Materials Science 

and 

G. B. McFadden and R. G. Rehm 
Mathematical Analysis Division 
Center for Applied Mathematics 

SUMMARY 

During solidification of a binary alloy at constant velocity vertically 
upwards, thermoaolutal convection can occur if the solute rejected at 
the crystal-melt interface decreases the density of the melt. We assume 
that the crystal-melt Interface remains planar and that the flow field is 
periodic in the horizontal direction. The time-dependent nonlinear 
differential equations for fluid flow, concentration, and temperature 
are solved numerically in two spatial dimensions for small Prandtl 
numbers and moderately large Sohmidt numbers. For slow solidification 
velocities, the thermal field has an important stabilizing Influence; 
near the onset of instability the flow is confined to the vicinity of the 
crystal -melt interface. For fixed velocity as the concentration increases, 
the horizontal wavelength decreases rapidly; a phenomenon also Indicated 
by linear stability analysis. The lateral lnhoaogenelty in solute concentra- 
tion due to convection is obtained from the calculations. For a narrow 
range of solutal Rayleigh numbers and wavelengths, the flow is periodic in 
time. 

Experimental observations of unidirectional solidification of succinonitrile 
containing ethanol have revealed an interaction between convective 
flow, Interface shape, and the localized development of a cellular Interface. 
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The effect is attributed to thermally inc ; ced convection and can be modified 
by changing the thermal field in the interface region. 

in collaboration with H. E. Gllcksman and Q. T. Fang of Rensselaer Polytechnic 
Institute, and R. F. Boisvert of the Scientific Computing Division, we have 
analyzed the stability of the parallel flow between a vertical crystal-melt 
Interface and a vertical wall held at a temperature above the melting point 
of the crystal. Three modes of instability occur: (1) a buoyant mode, 

(2) a shear mode, and (3) a coupled crystal-melt mode. For Prandtl numbers 
greater than approximately two, the coupled crystal-melt mode occurs at a 
lower Grashof number than the other two modes. These calculations and similar 
calculations for a cylindrical geometry were motivated by and are in general 
agreement with recent experiments on succinonitrile at Rensselaer Polytechnic 
Institute. 
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Introduction 


The general aim of this task Is the theoretical and experimental study 
of the fluid flow and solute segregation which occur during solidification. 
The nature of the fluid flow, its effect on the shape of the crystal-melt 
interface and on the resulting distribution of solutes is examined. In 
particular, the role of solutal convection (gravity dri . . low due to solute 
gradients) during the vertical directional solidification u binary alloy 
is considered both theoretically and experimentally. 

We present a disscussion of the numerical calculations first and then a 
description of the experiments: each Section can be read independently. 

The results of a number of calculations of the fluid flow and concentration 
field during directional solidification of an alloy with Prandtl number of 
0.01 and Schmidt number of 10 are described in the Numerical Results section. 
The experimental section describes the fluid flow and crystal-melt inter- 
face shapes observed during the directional solidification of succinonitrile 
containing small additions of ethanol. 

In an Appendix, we describe a linear stability analysis relevant to 
experimental work being conducted at Rensselaer Poytechnic Institute. As 
described in the Appendix an interface instability occurs under conditions 
for which in the absence of fluid flow, the crystal-melt interface would be 
morphologically stable. 
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NUMERICAL RESULTS 


The basic equations and methods for calculating the fluid velocity, 
solute, and temperature fields during dire- Clonal solidificati >n of a binary 
alloy vertically upwards have previously been reported [1]. In this report 
we present the results of a number of numerical computations for an alloy 
of Prandtl number P * 10“^ and Schmidt number Sc » 10. In general we will 
use dimensionless variables; Tables I and II define dimensional and dimension- 
less quantities, respectively. 

The solutal Rayleigh number Rs* at the onset of convective instability 
can be determined by linear stability calculations. For convective modes of 
instability, the results of the linear theory are approximately the same 
whether the planar crystal-melt interface is allowed to deform or Is maintained 
rigid. The critical solutal Rayleigh number as a function of V/D is plotted In 
Fig. 1. In these calculations we have kept the Imposed temperature gradient 

i 

in the liquid constant so that Ra(V/D)^ is constant; thus Ra decreases rapidly 
as (V/D) increases. Previous calculations [2-4] (see e.g. Table 3 of 
ref. 2) have Indicated that as V/D increases the thermal field does not affect 
tht critical solutal Rayleigh number. In fact, Hurle, Jakeman, and Wheeler [4] 
showed that for sufficiently large V/D the dependence of density on temperature 
could be completely neglected, i.e., they set the Ra - 0, and obtained results 
in agreement with those #f ref. 2, which corresponded to a liquid temperature 
gradient of 200K/cm. Thus, linear stability theory indicates that for 
small values of V/D, the thermal field is Important while for large values 
of V/D the thermal field is unimportant. 

The different behavior in these two regimes is further illustrated in 
Figs. 2 and 3. The critical bulk solute concentration C for the onset of 
instability is plotted as a function of the horizontal wavelength X for 
D/V ■ 0.015 and 0.15 cm In Figs. 2 and 3, respectively. The concentration 
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has been normalized by the minimim concentration C* as a function of wavelength; 
thus C/C* » Rs/Rs*. Linear stability theory predicts that the amplitude of 
a perturbation varies in time as exp [o r (X)t]. For a given value of 
C/C* > 1, there is a positive maximum value of a r as a function of X- The 
wavelength corresponding to the maximum value of a r for a given value of C/C* is 
indicated by the dashed line in Figs. 2 and 3. For D/V - 0.015, the wavelength 
corresponding to the maximum growth rate, slowly decreases as C/C* increases. 

In contrast, for D/V a 0.15 cm, the wavelength corresponding to the maxi mu a 
growth rate rapidly decreases as C/C* increases. 

He note that the unperturbed melt density decreases with height if 
Rs < kScRa/(P(k-l) ) or Rs < 429 Ra for the values given in Table II. This 
inequality holds for all of the non-linear numerical calculations presented 
in this report. Ho6t of our numerical calculations have been carried out for 
D/V “ 0.15 cm, where linear theory indicates that the thermal field is import- 

i 

ant and that the wavelength decreases rapidly with Increasing solutal Rayleigh 
number. For Rs/Rs*Hl we choose the horizontal wavelength as that 
corresponding to the onset of Instability as determined by linear theory. 

As the solutal Rayleigh number increases, we find that even though we start 
with a given wavelength, the actual numerical solution has a smaller wave- 
length, usually one half the original wavelength. We then use this new 
wavelength as the starting wavelength for calculations at larger solutal 
Rayleigh numbers. This behavior Is illustrated In Fig. 4. For D/V ■ 0.15 cm, 
the wavelength at the onset of convection X 0 ** 8.06 D/V. At Rs/Rs* 2 1.62 
and 1.65, the wavelength decreases to Xq/ 2 and X 0 /3, respectively. At 
Rs/Rs* 2 3.5, the wavelength decreases again from X 0 /3 to X 0 /6. As will 
be discussed, the numerical solution may depend on the initial conditions. 

Over a narrow range Rs/R we have found time dependent solutions which 
are indicated by the dashed line in Fig. 4. 


23 



Some testa of the numerical code have previously been described [1]. As a 
further test, we carried out calculations at \ - Xq and Rs/Rs* a 0.977 and 
i.020 on a (12, 36) grid and found no flow and flow, respectively, in agreement 
with linear stability analysis. 

Figures 5-9 illustrate numerical results for Rs/Rs* “ 1.29. In our numerical 
calculations, we use a dimensionless time, - tv/H^. We also define a dimen- 
sionless time T c based on the crystallization velocity, viz., T c « V^t/D, and 
note that T^/T^ 111 (VH/D)^/Sc. In the absence of convection, the time for the 
exponential concentration profile in front of the crystal-melt interface to 
develop is of the order of T c /k for k«l. Specifically, for k - 0.3 the 
Interface concentration reaches 95Z of its final value in a time corresponding 
to T c = 10.2 or = 5.9 for (VH/D) - 4.17 as in Figs. 5-9. For Rs/Rs* - 1.29, the 
time to approach steady state in the flow, solute, and temperature fields corres- 
ponds to = 14 as shown IniFig. 5 In which the maximum value of the stream 
function Is plotted against time. In Figs. 6-7 we show the steady state stream 
function and concentration respectively as a function of position. At the 
top boundary, we have assumed that the fluid velocity and perturbed temperature 
gradient vanish, and have used a mixed boundary condition on the solute 
concentration, which ensures solute conservation [5]. The top boundary is an 
artificial boundary, and we would like to place it sufficiently far from the 
Interface so that it does not influence the flow and solute segregation near the 
interface. From Fig. 6-7 It is evident that the top boundary is having little 
effect on the flow and solute concentration near the Interface. As the solutal 
Rayleigh number Increases, the flow is less confined to the Interface and even- 
tually it is not feasible In terms of computational times to move the upper 
boundary far enough from the interface. As illustrated in Fig. 8, the solute 
concentration *t the crystal -melt Interface varies by 6QZ. The average solute 
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concentration at the interface deviates by about 3Z from the correct value 
due to the finite mesh (9 x 36) used in the calculation. In Fig. 9 we show 
the average solute concentration in the liquid as a function of distance from 
the crystal-melt interface. For comparison, the exponential concentration pro- 
file which would hold in the absence of convection is also shown in Fig. 9. 

Over a narrow range of Rs/Rs* and X - Xq/ 3 , we have found 
time-dependent solutions. These are shown in Figs. 10-21; the calculations 
were carried out on a (12 x 72) mesh with (VH/D) “ 4.17 and Rs/Rs* - 3.3. 

Figs. 10 and 11 show the maxi mum value of the stream function and the 
average interface cone ntratlon in the melt, respectively, as a function 
of time. The period of oscillation is approximately ATjj 2 0.56 
or AT C s 0.97, i.e., the interface has moved a distance 0.97 (D/V) 
during one cycle. The magnitude of the average interface composition 
oscillation is small, viz. about 1.6Z, which Is considerably less than 
the transverse segregation of about 18Z. The magnitude of the stream 
function oscillation is considerably greater. Figures 12-21 show the 
stream function and concentration as a function of position at various 
times during one period of the oscillation. It is clear from these 
figures that the second vertical cell is undergoing the largest oscillation 
with the stream function varying in magnitude from 0.119 to 0.234. Bote 
that the flow does not reverse in direction; rather these are oscillations 
superimposed on a base flow. This may explain why linear stability 
calculations, assuming a motionless base state, have not revealed any 
time dependent phenomena in this parameter range. We have also observed 
that the steady state solution may depend on the Initial conditions, 

I.e., there may be multiple steady states. As Initial conditions, we 
have (1) perturbed a uniform concentration field, (2) perturbed a steady 
state (exponential) concentration field, and (3) used a solution for a 


25 





particular value of R a as the initial condition for a different value of Rs. 

In Figs. 22-23, we show the steady state stream function and concentration 
field, respectively, for Rs/Rs* ** 3.6 using the oscillatory solution at 
Rs/Rs* - 3.3 (see Figs. 10-21) as the Initial condition. In contrast. 

Figs. 24-25 show the steady state stream function and concentration 
field, respectively, for exactly the same conditions as Figs. 22-23 
except that the initial conditions consisted of a perturbed exponential 
concentration field. For this case, the final solution has one half the 
wavelength of the initial condition. 

In Figs. 26-2 7 , the stream function and concentration field, 
respectively, for Rs/Rs* « 5.03 and (VR/D) - 4.17 are shewn. Figs. 28-29 
are identical to Figs. 26-27, respectively, except that (VH/D) and the 
number of mesh points in the vertical direction have been doubled. It is 
evident in Figs. 26 and 27 that the effect of the upper boundary is not 

i 

negligible. As can be seen In Fig. 28 with (VH/D) - 8.33, the stream % 

function vanishes at a distance of approximately 5.4 (D/V) from the 
interlace. The solute segregation at the interface is about the same in 
Fig. 27 and 29, viz., 15Z. 

Figures 3-29 were for D/V - 0.15 cm; in Figs. 30-33 we show the results 
of a calculation for D/V - ff.015 cm, Rs/Rs* - 1.325 and Ra - 0.127, corres- 
ponding to the same temperature gradient a a in Figs. 3-29, As previously 
discussed, we do not expect the temperature field to stabilize the flow 
since the Rayleigh number is very small. As shown in Figs. 30-31, even 
though we have taken H * 20 (D/V), the flow is influenced by the upper 
boundary, and is certainly not confined to the vicinity of the crystal-melt 
interface. There are clearly fairly large lateral solute gradients far 
from the interface. Fig. 32 shows the solute concentration at the crystal-melt 
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interface; the solute segregation is about 60Z. The average solute 
concentration as a function of distance froa the interface is compared 
with the solute profile In the absence of convection in Pig. 33. The 
average Interface concentration and the average concentration gradient at 
the interface are determined by the boundary couditions, and are i^iipevsut 
of the fluid flow as can be se~n from Pig- 33- However, at about 2 D/V 
froa the interface, the average solute concentration with convection is 
about 40Z higher than the unperturbed (convectionless) concentration. 

Since the Rayleigh number was mail for the results of Pigs, 30-33, we 
carried out an Identical calculation with Ra - 0 to determine whether the 
thermal field can be completely neglected. The maxima value of the 
streaa function was 80 Z larger for Ra - 0 than for Ra » 0,127 indicating 
that the thermal field plays a role even for Ra as small as 0,127. 

The numerical results presented here provide Insight Into the solutal 

* - 

convection which can occur during directional solidification. The range 
of parameters which can be studied is severely limited by the computational 
times required for these calculations. Results for different values of 
the distribution coefficient and temperature gradient in the liquid would 
be desirable. Ve have limited the calculations to s Schmidt number of 
ten, which is appropriate for semiconductors . The time to reach steady 
state appears to be proportional to the Schmidt number for large Schmidt 
numbers so that calculations for metallic melts with Schmidt numbers of 
the order of 100 will take ten times as much computer time. The interaction 
of the solutal convection, which is driven by the adverse concentration 
gradient near the crystal-melt interface, with the thermal convection due 
to horizontal temperature gradients in a real furnace is an important 
area for further research. 
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EXPERIMENTAL 


Background 

Succinonitrile containing ethanol provides a useful system for the 
investigation of thenao-solutal effects during solidlf ication. The 
density of ethanol being significantly less than that of succinonitrile, 
unidirectional upward solidification in this system leads to solute 
distributions which are predicted to be unstable with respect to convective 
perturbations- We previously determined the succinonitrile-rich end of 
the phase diagram for this system and, on the basis of the measured phase 
diagram parameters and the published physical properties of succinonitrile, 
predicted the critical concentration of ethanol in succinonitrile above 
which instability occurs, as a function of the velocity of upward 
solidification- The predicted instabilities are either interfacial, in 
which case the wavelength Is short and the shape of the solid-liquid 
interface deviates from planarity, or convective, in which case the 
wavelength is longer and the liquid above the interface starts to flow. 

At a fixed temperature gradient, the concent ration of ethanol required to 
produce interfacial Instability decreases with increasing solidlf ication 
velocity, while the concentration required to produce convective instability 
increases with increasing solidification velocity. At some compositions, 
there therefore exists a range of solidification velocities in which 
unidirectional upward solidification is predicted to he stable; above 
this velocity range lnterfaclal Instability is predicted and below this 
velocity range convective instability is predicted. 

Our experimental observations of solidification processes in the 
region of these transitions have used both pure succinonitrile and 
succinonitrile doped with an appropriate concentration of ethanol to 
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investigate the behavior when a finite range of stable velocities Is 
predicted. Because of the low solute distribution coefficient for ethanol 
in succiaonitrile (k ■- 0.044), concent rations of interest are very low 
(2 - 5 x 10*“^ wt.X ethanol) and the samples were thus made from ultra- 
purified succiaonitrile and sealed under vacuum. 

Perfectly uniform upward solidification of a substance with a normal 
thermal expansion coefficient (I.e., one which expands with Increasing 
temperature) is convectively stable if no solutes are present. The 
convective instability predicted for the succinonltrile-ethanol system 
arises from the rejection of the lower-density solute by the advancing 
interface. In a real experiment , however, there are sources of convection 
which may be described as purely thermal and which Interfere with 
observation of the effects predicted for truly unidirectional solidification. 
Any factor which causes the lpotherms in the liquid to deviate from 
horizontal planes will be a driving force for convective flow. 

There are two important causes for thermally Induced fluid flow in 
the liquid above a solidifying Interface when unidirectional solidification 
of a material in a tubular container is attempted by drawing the container 
downward through a temperature gradient. First, heat losses or gains 
through the container walls to the furnace or outside environment 
necessarily imply the presence of a radial component in the temperature 
gradient. These mu at be minimized by design of the crystal growth 
apparatus and by control of the Interface position. The second origin 
for thermally Induced convective flow is tbe distortion of the Isotherms 
by the container walls. When no solidification is taking place, there is 
a gradient discontinuity at the solid-liquid Interface due to tbe different 
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thermal conductivity of the solid and liquid. Because this gradient 
discontinuity is not present in the container walls, the isotherms in the 
region where the interfce meets the walls must be distorted from planarity. 
Succinonitrile is unusual in that the thermal conductivity of solid and 
liquid are, within experimental uncertainty, equal, so that at zero 
solidification velocity this effect is not important. When solidification 
starts, however, a gradient discontinuity Is produced at the solid-liquid 
interface by the emission of latent heat of fusion. The heat emitted per 
unit area of the interface is vL where v is the velocity and L is the 
latent heat (4.7 x 10^ J/zo} for succinonitrile). Because in a typical 
experiment the discontinuity produced by this latent heat is not negligible 
compared to the externally imposed gradients, a significant thermal field 
distortion is produced by the container walls. 

For the unidirectional solidification of metals it is often possible 
to use containers having ouch lower thermal conductivity than the metal, 
in which case the thermal distortions caused by the container walls will 
be small. However, lateral heat losses to the environment can be large in 
metallic systems due to the high melting temperatures. For an organic 
material such as succinonitrile, however, the thermal conductivity of the 
container cannot ha ouch lower than that of the material and container 
wall effects are unavoidable. Radial heat losses to the environment, 
however, can be small because of the low melting temperature. Thus 
distortion of the Isotherms may be similar in metallic and organic systems, 
although the most Important cause of the distortion may differ. 

In our experiments we have Investigated the effects of radial 
gradients on “unidirectional" upward solidification under conditions 
close to thoso predicted to result in solute-induced Interfacial or 


30 


convective instability, and have investigated methods to reduce the 
effects of radial gradients. We have also measured densities of 
succinonitrile-ethanol mixtures because earlier calculations were based 
on estimated values. 

Experimental Methods 

a. Density Measurements 

The density of pure liquid succinonitrile and succinonitrile containing 
four different concentrations of ethanol (up to 6.67 wt.Z) was measured 
by determining the temperatures at which calibrated standards with 
densities of 0.9700, 0.9800, and 0.9850 had neutral bouyancy. Distilled 
succinonitrile was used and measurements were carried out under air at 
atmospheric pressure. The temperature of neutral bouyancy was determined 
to within 0.i°C. 

b. Convective Flow During Unidirectional Solidification 

Pure succinonitrile and succlnonltrlle-ethanol samples are contained 
within borosllicate glass tubes 19 mm in diameter and approximately 50 ca 
long. The residual impurity content of the succinonitrile samples, before 
doping with ethanol, va a found to be equivalent to approximately 10"^ 
wt.Z (lppm) of ethanol. In terms of its effect on the solidus and llquldus 
temperatures. Most experiments were carried out on a sample containing 
2.6 x 10~3 wt.Z ethanol. 

To produce solidification, the sample tubes are drawn downward through 
a temperature gradient established by a thermostatically controlled lover 
cold water jacket and an upp'T heater which is either a similar hot water 
jacket or an electric heater. Growth velocities can be from 0.4 to 10 ixn/s. 



An important constraint for these experiments is the requirement of 
good visibility in the interface region. This can be attained only if 
the cylindrical tube is surrounded in the interface region by a medium 
with a refractive index similar to that of the succinonirrile, but which 
itself has a flat outer wall for microscopic observation. The best image 
is obtained if the index-matching medium is a fluid contained in a flat- 
walled outer container, but such a system could also undergo convective 
processes which would tend to perturb the temperature field within the 
solidification tube. 

Ve have therefore surrounded the interface region of the solidification 
tube with close-fitting transparent solid pieces eith flat outer windows. 
Heaters have been incorporated into these pieces in the region of the 
solid-liquid interface to provide control over the isotherm shape within 
the sample tube. The heaters are in the form of a flat (pancake-like) 
coil so that they can duplicate in the transparent solid the temperature 
discontinuity existing in the sample at the solid-liquid interface. The 
required heater power is proportional to the velocity of motion of the 
sample: for a velocity of 1 ym/s, a heater power of about 0.5 ¥ is 

appropriate. 

The sample tube is photographed from the side by a horizontally mounted 
microscope. To measure convective flow, small neutrally bouyant latex 
particles are incorporated into the sample. These are photographed using 
a triple flash exposure with unequal intervals between the flashes, to 
allow determination of the direction of particle motion. The particles 
are seen only in dark field illumination, but the Interface shaps is seen 
better in bright field illumination. 



Results 


a. Density 

The aeasured densities of succinonltrile-ethanol solutions are shovn 
In Figure 34. The aeasured densities are close to those which are calculated 
under the assumption that no change of total volume takes place when the 
two coaponents are mixed, although a deviation froa this behavior alght 
have been expected on the basis oi the liquid-phase miscibility gap In 
the succinonltrile-ethanol system. The lines of constant density lie 
almost exactly parallel to the liquldus line In the composition and 
temperature range which was measured. The density of the liquid at the 
liquldus is thus constant, 0.9880 g/ca^. This leads to the unusual result 
that the liquid at the surface of a sucdnonltrlle crystal growing Into a 
succinonltrile-ethanol mixture has a constant density, no matter what the 
shape of the crystal, so long as the ethanol concentration does not exceed 
approximately 7 wt.Z. Liquid not Immediately at the surface, however, is 
not constrained to be at the liquldus and therefore has a variable density. 

The measurements Indicate a concentration dependence of liquid 
density a c •* -<3p/3c)/p ■ 2.7 x 10 - ^ (wt.X)"^ which is close to the 
value used previously (3.07 x 10”^ (wt.X) - ^. 

b. Convective Flow During Unidirectional Solidification 

When the sucdnonltrlle sample is held stationary In the directional 
solidification apparatus, the solid-liquid Interface is horizontal and 
very slightly concave upwards. The amount of this concavity, which 
probably results froa a snail radial heat loss, depends upon the position 
of the Interface between the hot and cold jackets, the nature of the 
surrounding medium, and other details of the growth apparatus design. A 



concave interface indicates that the liquid in the center of the tube is 
warmer than that near the wall at the same level. Therefore, as expected, 
there is an upward flow of liquid in the center and a downward flow closer 
to the walls, with the liquid immediately above the solid-liquid interface 
moving radially inward. The upward flow velocity at the tube center is 
typically 10 jjm/s. The only microstructural features present at the 
interface in either pure or ethanol-doped samples are faint grooves due 
to grain boundaries or sub-boundaries. 

When solidification is produced by drawing the sample tube downward, 
the solid-liquid interface becomes strongly concave in either pure or 
ethanol-doped succinonit rile. Simultaneously there is a substantial 
increase in the velocity of upward flow at the tube center (to about 
40 \mJ s for a growth velocity of 2 ^n/s). The pattern of the flow 
remains unchanged. 

These flow velocities are sufficiently large that they can be expected 
to seriously distort the solute distribution in the liquid above the 
solid-liquid interface. It was found, however, that the pancake-shaped 
auxiliary heater surrounding the interface region could be effective at 
altering this flow pattern. 

When the auxiliary heater is not used, the convective flow can result 
in a transient macroscopic interface shape change which then leads to a 
localization of interfacial instabilities. The sequence of events which 
constitute this phenomenon is as follows: 

a. The interface is equilibrated in its almost planar configuration 


at zero growth velocity. 



b. Growth is started at a velocity someV'-t less than that which is 
expected to produce interfacial instabilities. 

c. The interface becomes concave upwards as a result of the heat 
flow effects mentioned above. 

d. Following an incubation time which depends on the growth velocity, 
an additional deeper depression develops in the center of the concave 
interface. The size and shape of this depression also depends on the 
growth velocity. 

e. The depression deepens rapidly and interfacial instabilities 
develop within it. 

f. The cellular interface within the depression advances rapidly 
and catches up to the general contour of the concave interface. 

g. The concave interface then propagates at steady state with a 
central cellular region in tfc$ area where the depression had been. 

The incubation time for development of the depression is longer at 
lower growth velocities. At 1 um/s, the depression has not been 
observed, and the interface remains bowl-shaped but smooth (no interfacial 
instabilities). At 2 ym/s the conditions for formation of the 
depression are apparently marginal and the incubation time in different 
runs was found to be between 26 and 125 minutes. At a growth velocity of 
3 vin/s the incubation time was 13 to 15 minutes. At 4 \m/a the 
incubation time becomes more difficult to define because of shape changes 
(see below) which wake the pit less distinct from the overall concavity 
of the Interface. 

The transient interface pit takes its most dramatic form under the 
most marginal conditions for its formation. Thus in the 2 ^m/s run in 
which the pit took 125 minutes to develop, the pit was very deep, narrow, 
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and pointed before a cellular structure developed near its bottom and 

spread upwards along the walls (Figure 35 a). At faster growth rates, or 

when the growth conditions at 2 u®/s were soqxwhat different and the 

pit formed more rapidly, the pit was broader, slialiower, and more rounded 

when the cellular structure developed within it (Figure 35 b). After the 

cellular interface formed within the pit, cells spread rapidly up the 

sides of the pit as far as the point where the pit merged into the overall 

bowl-shaped interface. The pit region then filled in rapidly with cellular 

material until the leading edge of this material was essentially coincident 

with the overall bowl-shaped contour of the interface. Fig' re 36 shows 

the development of the pit in a run at 2 p*./s. The position, relative 

to a fixed point in the apparatus, of the bottom of the pit is shown; 

after the cellular substructure develops, the position o. the first solid 

(the leading edge of the cells) and the last liquid (the last entrained 

» 

intercellular liquid) are shown. 

After the development of the cellular interface in the central pit, 
the convective flow pattern in the liquid above this region is different. 
Many particles are seen indicatiag fluid flow downward into the central 
cellular interface. The flow appears to be less regular than the simple 
pattern which prevails before the cells develop, but we have not yet 
discerned its detailed structure. x 

When the auxi liar>*iieat-er Surrounding the interface is- activated, 
the formation of the central pit can be suppressed. Ideally, the heater 
would be located exactly at the interface level and would transfer heat 
uniformly to the sample container walls from all sides. In practice it is_ 
difficult to arrange that the interface lies exactly in the plane of the 
heater because the vertical position of the interface is affected by many 
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factors including the growth velocity and the auxiliary heater power. 
Circumferentially uniform transfer of heat to the container walls can 
also only be approximated. Nonetheless, auxiliary heaters are very 
effective at changing the convective flow pattern and eliminating the 
central pit. 

As an example, on one run at 2 ixt/s a central pit developed 33 

minutes after the start of growth and the Interface within the pi then 

became cellular. The auxiliary heater was then turned on and its power 

level adjusted to give, as closely as possible, a flat Interface near the 

walls of the sample tube. After approximately 30 minutes, the central 

cellular region of the interface had disappeared. The interface at this 

point was relatively flat but tilted to one side: the convective flow in 

the liquid was downward on one side of the tube (where the tilted Interface 

was higher), across the tube above the interface, and upward on the 

» . 

opposite side (where the Interface was lower). Similarly, if growth is 
started when the heater is activated to an appropriate level, the central 
pit does not develop. 

Discussion 

Thes^* experiments have demonstrated the importance of a previously 
unrecoreded interaction between convective flow, interface shape, and the 
development of cellular interfaces. As observed in the succinouitrlle- 
ethanol system, the effect occurs over a small range of growth velocities, 
but this range of growth velocities is an important one because it 
represents the conditions which might be chosen for an attempt to grow a 
crystal of homogeneous, cell-free material. The predicted velocity for 
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interfacial instability in the system used for these experiments is 
approximately 10 ym/s (depending on the value of the temperature 
gradient, which is close to 10 K/ cm), and the interactive effects were 
observed at growth velocities of 2 to 4 ;m/s. In addition to being 
limited to a narrow velocity range, the effect is transient in nature, 
with the result that the evolution cf interface shapes by this process 
will not be understood if observations are not sufficiently frequent. 

Figure 37 summarizes the nature of these interactions in pure and 
solute-doped systems. At zero growth velocity, the solid-liquid interlace 
in either system can be planar and horizontal. Vhen solidification 
occurs, the thermal effect of the container walls results in a concave 
solid-liquid interface and a convective flow pattern in which the liquid 
above the solid-liquid interface sweeps radially inward and then up along 
the tube center. 

In the pure material, the interface remains smooth and bowl-shaped* 

If solute is present, the radial inward flow redistributes the solute 
rejected by the interface and the concentration at the interface is 
r^ghest in the central region. The resultant local depression of the 
melting point in this region causes the formation of the pit, which 
becomes progressively deeper as the geometry for diffusive and convective 
removal of solute becomes less favorable. 

The absence of the pit in the pure system, and the fact that the pit 
formation can be suppressed by the auxiliary heaters, both support the 
interpretation that the pit formation is the result of lateral redistribution 
of solute at the solid-liquid interface. Formation of a macroscopic pit 
is possible only if the interface does not immediately become morphologically 
unstable, and the observed velocity dependence meets this condition. 
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iable 1. Definition of syiab used for physical quantities. 

v (kinematic viscosity) 
k (thermal diffuaivity of melt) 

k b (thermal diffuaivity of crystal) 

D (diffusion coefficient of melt) 

g (acceleration of gravity) 

a (thermal coefficient of expansion of melt) 

(Xq (solutal coefficient of expansion of melt) 

k L (thermal conductivity of melt) 

kg (thermal conductivity of crystal) 

k (distribution coefficient) 

Co. (bulk concentration of solute) 

C* (bulk concentration corresponding to onset of convection) 
(temperature gradient in melt) 

H (melt height) 

V (crystal growth velocity) 

X (horizontal wavelength) 
t (time) 


* 


1 


\ 
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Table 2. Dimensionless quantities used In calculations (unless otherwise 
specified) 


Ratio of Thermal Conductivities 

Ratio of Thermal Dif fusivities 

Distribution Coefficient 

Schmidt number 

Prandtl number 

Thermal Rayleigh number 

Critical solutal Rayleigh number 

Dimensionless Crystal Growth Velocity 

Ratio of solutal Rayleigh number, 

to critical solutal Rayleigh number 

Ratio of melt height to horizontal wavelength 


k L /k s -1.0 

k/cj - 1.0 
It - 0.3 

Sc - V/D - 10.0 

P - V/k - 0.01 

8a - (g<iG L /VO(D/V) 4 - 1268. 

r£ - (go c cl/VD)(D/V) 3 - 963. 

VH/D - 4 - 20 

Rs - (8a c <VVD)(D/V) 3 , 

8s/8S - VARIOUS VALUES 
HA - VARIOUS VALUES 
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RAYLEIGH NUMBER 


The solutal Rayleigh number, Rs, at the onset of convection as a 
function of the ratio of interface velocity, V, and the diffusion 
constant, D. The temperature gradient in the melt was kept 
constant so that the thermal Rayleigh number, Ra , varies as the 
inverse fourth power of V/D. The values of the other parameters 
are given in Table 2. 





c/c* 


LINEAR THEORY FOR D/V - .015 CM 


o 



Pig. 2 The normalized critical bulk solute concentration, C v at the onset 
of convection as a function of the horizontal wavelength (in units 
of D/V ■ 0*015 cm.) of a sinusoidal perturbation* The wavelength 
corresponding to the waxirsua growth rate of the perturbation is 
indicated by the dashed line. 
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Fig* 3. The normalized critical bulk solute concentration, C, at the onset 
of convection as a function of the horizontal wavelength (in units 
of D/V * 0.15 cm.) of a sinusoidal perturbation. The wavelength 
corresponding to the maximum growth rate of the perturbation is 
indicated by the dashed line. 
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NONLINEAR CALCULATIONS FOR D/V 
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Fig, 4. The vertical lines indicate the wavelengths and solute concentrations 
for which nueer leal calculations were carried out. The horizontal 
Hoes Indicate values of solute concentrations for which the wave- 
length decreases. Ttie solid curve is the same as in Fig, 3, 
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STREflM FUNCTI0N CHNT8UR MAP AT T = 24.32 
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Fig. 6. The steady state stream function as a function of position for 
the same parameters as in Fig. 5. Each vertical and horizontal 
division indicates the mesh used in the numerical calculations; 
although the calculation is done for one half wavelength in the 
horizontal direction, slightly more than one wavelength is show 
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CONCENTRATION 




Z 

Fig. 9. The average steady state concentration in the melt as a function 
of the distance (in units of D/V) from the crystal-meJt interface 
for the same parameters as in Fig. 5. The dashed line indicates 
the concentration in the absence of convection. 
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STREAM FUNCTION 




Fig* 10. The maximum of che stream function as a function of time (in units 

of H^/v) for Rs/Rs =» 3.3, VH/D D A. 17, and H/X « 1.551. The 4 

values of the other parameters are given in Table 2. The stream 
function is periodic in time; the spatial variation of the stream i 

function during one cycle is shown in Figs. 12-16. j 

i 

4 
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The average concentration at the crystal-melt interface as a function 
of time for the same parameters as in Fig. 10. The spatial variation 
of the concentration during one cycle is shown in Figs. 17-21. 
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CSNCEHTRflTJW C2NTBJR MflP AT T = 18.93 



Fig. 19. The solute concentration as a function of position at dimension- 
less time 16.93 for the same parameters as in Fig. 10. 
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The steady state * " 
Rs/Rs - 3.6, VH/D = 
The Initial state f 
concentration field 
wavelength of the i 
that the steady sta 


n function as a function of position for 
17, and H/X « 1.552 (same as Figs. 22-23). 
this solution was a perturbed exponential 
The steady state solution has one half the 
ial state. Comparison with Fig. 22 indica 
solution depends on the initial conditions 
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Fig. 25. The steady state concentration as a function of position for 
the same parameters as in Fig. 24. 


j 


66 











Pig. 28. The steady state stream function as a function of position for 
Rb/r 4 - 5.03, VH/D = 8.33, and H/X - 6.208. The values of 
the other parameters are given in Table 2. This calculation 
is identical to that of Fig. 26 except that the height has 
been doubled. 
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Pig, 29. The steady state concentration as a function of position for the 
same parameters as in Pig. 28. 
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Fig. 33. The average steady state concentration in the melt as a function 
of the distance (in units of D/V) from the crystal-melt interface 
Tor the same parameters as in Fig. 30. The dashed line 
indicates the concentration in the absence of convection. 
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100 120 140 160 »{S0: 

TIME (minutes) 

Fig. 36. Interface positions during the formation of a '^ntral pit, c.s a 

function of time following the start of growth, 2 yra/s. The position 
of the bottom of tha pit is plotted until the cellular structure develops 
the leading and trailing ed;;£s of the cellular structure are then shown. 








APPENDIX 


Coupled Convective Instabilities at Crystal-Melt Interfaces 

S. R. Coriell, G. B. McPadden, and R. F. Boisvert 
National Bureau of Standards 
Washington, D.C. 20234 

and 

M. E. Glicksman and Q. T. Fang 
Materials Engineering Department 
Rensselaer Polytechnic Institute 
Troy, NY 12181 

The interaction between fluid flow and the crystal-melt interface is 
of central importance in determining the solute distribution in the 
solidified material, which in turn determines its properties* The effect 
of convective flow on morphological stability has been recently reviewed 
£ 1 J - Calculations by Delves [ 2 1 shoved that a flow parallel to the 
crystal-melt interface could increase stability and give rise to traveling 
waves on the interface* The Interaction of morphological instability 
with a aelt subject to thermosoJutal Instability has also been studied 
[3J • In both these cases, morphological instability occurs even in the 
absence of fluid flow, although in general the critical value of the 
parameter for the onset of instability is changed by the flow. 

Recent experiments [4-7 J have demonstrated an interface Instability 
un, _r conditions for which in the absence of flow, the crystal- aielt 
interface would be morphologically stable. In these experiments, a long 
vertical cylindrical sample of high purity succioonltrile was heated by 
passing an electrical current through a long, coaxial, vertical wire, so 
that a vertical melt annulus formed between the coaxial heater and the 
surrounding crystal-melt interface. The outer radius of the crystal was 
maintained at a constant temperature below the melting point (56«1°C). 
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With this arrangement the temperature increases in the melt and decreases 
in the crystal with distance from the crystal-melt interface, and 
consequently the interface would be morphologically stable in the absence 
of fluid flow. 

Buoyancy dictat- a that the f^uid flows upward near the heated wire 
and downward near the crystal-melt interface. There are analytic solutions 
of the fluid flow and temperature equations in which the flow velocity is 
vertical and the flow velocity and temperature are functions of the radial 
coordinate alone. Choi and Korpela [8] and Shaaban and Ozlslk [9J have 
calculated the critical Grashof number for axisymmetric instabilities of 
the axisymmetric flow between two vertical, infinite, coaxial, rigid 
cylinders held at different temperatures. We [6] have repeated these 
calculations by a different numerical technique and for physical properties 
appropriate to succinonitrile (Prandtl number of 22.8). The flow betweeu 
two vertical infinite coaxial cylinders containing succinonitrile is 
unstable to an axisymmetric perturbation above a Grashof number of 2150 
and the wave speed of this perturbation is comparable to the unperturbed 
flow velocity. 

When these calculations are carried out with the outer rigid cylinder 
replaced by a crystal -me It interface, quite different results are obtained. 
The onset of instability occurs by an asymmetric mode at a Grashof number 
of 176, which is in agreement with experiment. The linear stability 
calculations indicate that the density change on solidification and the 
crys tai^melt surface tension are of little Importance in determining the 
instability. The simulations predict that the period of oscillation 
(time for Che wave to traverse one wavelength) is proportional to the 
fifth power of Che gap width (the difference between the radius of the 
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crystal-melt interface and the wire radius). This has been verified 
experimentally with the period of oscillation ranging from 300 to 60000 8. 

In order to understand better the “coupled mode" of instability, we 
use linear stability theory to analyze a simpler geometry, viz., we 
consider the stability of the parallel flow between a vertical crystal- 
melt interface and a vertical wall held at a temperature above the melting 
point of the crystal. Three modes of instability occur: (1) a buoyant 

mode, (2) a shear mode, and (3) a coupled crystal-melt mode. The first 
two modes are well known from previous research [10] on the instability of 
the flow between two vertical walls. The behavior of these modes as the 
Prandtl number varies will also be described. 

Theory 

We consider an (x, y, z) Cartesian coordinate system such that in the 
unperturbed state the region x < 0 is crystalline, the region 0 < x < L 
is molten, and at x - L tb^re is a rigid isothermal wall, whose temperature 
is greater than the melting point of the planar crystal by an amount AT. 

We assume that the force of gravity acts in the negative y direction and 
that all quantities are independent of the z -coordinate. We measure 
lengths in units of L; temperature T in units of AT; time t in units 
of L^/v, where v is the kinematic viscosity, and velocities in units of a 
reference velocity D, which will be specified subsequently. The 
dimensionless equations for the fluid velocity u and the temperature are, 
in the Bous8ineaq approximation 

V- u - 0, (1) 

( I/3t) + Re (u*V)u - - Vp + V 2 ^ + (G/Re)Tj, (2) 

(3T/3t) + Re u-VT - (V 2 T)/P, (3) 
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where Re a UL/v is the Reynolds number, p is a dimension J ec.s , 

G - goATL-tyv^ is the Grashof number, j is a unit vector in the y- 
direction, P is the Prandtl number, g is the gravitational acceleration, 
and a is the coefficient of thermal expansion. The temperature Tg in 
the crystal satisfies 

(3T s /3t) « V^r/Pg, (4) 

where Pg «=■ ?</ eg « v/icg ai *d < aad <s are t ^ xe thermal dif fueivlties of 
the melt and crystal, respectively. 

The boundary conditions at the crystal-melt interface are that the 
tangential component of the fluid velocity vanishes and that the normal 
component satisfies 

u *n a - e(vn) . (5) 

Here n is the normal to the crystal-melt interface, e ® -1 + Ps/PL» 

PS and pl are the crystal and melt densities, respectively, find v 
is the local solidification velocity. At the crystal-melt interface, the 
temperature is continuous and equal to the equilibrium temperature , viz,., 
T S I - T LI - -TK, (6) 

where the sul 3cript 1 denotes evaluation at the crystal-melt interface. 
Temperature is c^asured relative to the dimensionless melting point T^/AT 
of the nlanar crystal-melt interface, T m T^ff/ (LH V AT) , y is the 
crystal-melt surface free energy, H v is the latent heat of fusion per 
unit volume of solid, and K is the mean curvature of the crystal-melt 
interface. Conservation of energy at the crystal-melt interface requires 

A(v*n) - (-k r VT + 7T a )«n, (7) 

where A « vH v Re/k 3 AT and k r is the ratio of the thermal conductivity 
of the melt k^ to that of the crystal kg. 
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At the rigid wall y “ 1, • 

u - 0 (8) 

and 

T - 1. (9) 

In the crystal we require that far from the crystal-melt interface, the 
teo^erature gradient is constant and assumes a value such that the 
unperturbed planar Interface is stationary. 

A solution of the differential equations satisfying the boundary 


conditions is 

Uyo - [1/4 - (x - l/2) 2 ][(G/6Re)(x-l/2) + 61], (10a) 
T 0 - x, (10b) 
Tso “ M. (10c) 
Po - [G/2Re - 12I]y, (lOd) 


where I is the net flow in the y-direction, i.e. , the integral of Uy from 
x - 0 to x - 1. The subscripts x aud y denote the x and y components of 
u and the subscript o denotes the unperturbed or base flow and tenperature 
fields. 

We examine the stability o'i the flow with respect to small 
perturbations. We decompose the fluid velocity, temperature fields, and 
crystal-melt Interface shape xj Into an unperturbed and perturbed part, 
and assume that the y and t dependences of the perturbed quantities are 
of the form F(y, t) * exp (iwy + at), where ai Is a spatial frequency and 
o ■ a r + ia^ is a complex number which determines the temporal 
behavior of the system. If a r > 0 for any value of u, the system is 
unstable. We write 
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“ W x (x)F, 

(11a) 

Uy ° UyQ + Wy(x)F f 

(lib) 

T - x + T(x)F, 

Ole) 

T s = k r x + T 8 (x)F, 

(lid) 

p° p 0 + pi(x)f. 

(lie) 

XI ■= 6F , 

(Ilf) 


where 6 is the amplitude of the perturbation of the crystal-melt Interface. 

The linear ordinary differential equations for the perturbed 
quantities, W x , Vy, pi, T, and T s can be written as 


DW X + iwVy » 0, (12a) 
\W X + ituDWy + Dpj ° 0, (12b) 
-Re(Du yo )W x + (D 2 - A}W y - iwpj + (G/Re)T = 0, (12c) 
-PReW x + (D 2 - A t )T “ 0, (12d) 
(d 2 _ u 2 _ P s o)T s = 0 , (12e) 


where D ° (3/3x), A ° w 2 + o + iuReu yo> and - u 2 + P(o + itoReu^). 

The perturbed temperature field in the crystal, can be found analytically 

and is of the form exp + PgcO^^x]. 

Tt is desirable for purposes of numerical computation to prescribe the 
boundary conditions at x - 0 rather than at the crystal-melt interface x « x^. 
This is easily accomplished in linear theory because any function, e.g., f(x), 
evaluated at x 3 xj - Ar can be written as f(xx) “ f(0) + 6F(3f/3x), where 
the partial derivative is evaluated at x - 0. Upon eliminating Tg and 6, 
we find the following boundary conditions for the fluid variables, W x , 

Wy, and 1 at x « 0: 

(1 + rw 2 )ReW x - eaT » 0, (13a) 

(1 + ru; 2 )W y - ( Du y o ) T « 0, (13b) 

[(u 2 + P c o) 1/2 (k + Tu) 2 ) + aA/Re] T - k (1 + I\d 2 )(3T/3x) - 0. (13c) 

c* r r 
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These boundary conditions are furhter simplified by taking T ■ c « 0. 

It is not surprising and numerical calculations have verified that these 
are excellent approximations. 

The boundary conditions at the rigid isothermal wall at x «* 1 are 

V x - Vy - T - 0. (14) 

In order to understand the physical mechanism of instability, it is 
useful to compute the average kinetic energy of the perturbed flow lield. 

The average kinetic energy in the disturbances is obtained by integrating 
the quantity (Real u^)^ + [ReaKu^ - Uy Q )] 2 , as given by eq. (11), over a 
basic cell having unit width in ihe redirection and one wavelength 
(2r/a) in the y-direction (Real indicates the real part of the complex 
functions u x and Uy) . At the onset of instability (o r - 0) the rate 
of change of the kinetic enorgy in the disturbance vanishes, and from 
cqs. (I2b) and (12c), it follows that 

B 

-<K | 2 + J^yj 2 + “ 2 | W x | 2 + “ 2 ] w yl 2 + (G/Re) (*««UW y *)> 

- Re(DUy 0 )(Real(V ; Wy*>> - Real(tf y (0)DW y *(0) ] - 0, (15> 

where the < > Indicates integration from 0 to 1 in the indirection and 
the * indicates the complex conjugate. Ihe sum of the iirst four terms 
repres nts energy loj< io viscous dissipation e r 1 1 j always negative. It 
Is balanced by tue &\ — of the remain ig three torms. Ve denote the terms 
proportional tr G and (Duy 0 ) as buoyant and shear terms, respectively. 

The remaining term VyDWy evaluated at x 0 is due to the crystal-melt 
interface arid would vanish for a rigid interface. In the results to be 
presented subsequently, we have found modes of instability for which one 
of the terms, i.e., buoyant, shear, or interface, is dominant. 
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The numerical methods used to solve the linear • : * > a -in 

defined by the differential equations leas* 12a-4] and the V. udniy conditions 
[eqs. 13-14] are all similar to those described previously [3]. We 
usually vary G and keeping a r and all other parameters constant, 
until a solution of the differential equations is found which satisfies 
uhe boundary conditions. The program SUPORT [II] is used to solve the 
differential equations, and the program SNSQE [12-13] is used for the 
non-linear iteration procedure. Calculations of the critical Grashof 
number for the flow between two rigid plates are iu agreement (to 4 
significant figures) with Table 3 of Korpela, Gozum, and Baxi [14] • 

Numerical Results 

For the calculations reported here, we have taken 1*0 (ste eq. 10a) 
so that the net flow in the y-direction vanishes, and we choose our 
reference velocity U such that Re n G. Further we take T a e D 0 

i 

and k r “ P/P 5 » 1, i.e. , we neglect both the crystal-melt surface energy 
and the density change during solidif ication, and we assume that crystal 
and melt have the same thermal properties. In mo 6 t of our calculations 
we set A « vRyRe/kgAT - H v gaL^/k s v ° 1.017(10^), which corresponds to 
succinoni t rile with L *» 0.25 cm. The corresponding value of A for 
lead with L = 0.25 cm is A * 625. In general we calculate the onset 
of instability viz. o r - 0. For the Prandtl number P and all other 
parameters fixed, we calculate the Grashof number C and as functions 
of the spatial frequency ui for o r 0 0. Typical results fos succinonitrile 
P « 22. S are shown in Figures 1-4. As previously discussed the instabilities 
can be characterized as buoyant, shear, or interface. For P a 22.8, the 
mode with the lowest Grashof number shown in Figure 1 is primarily an 
interface mode. The next highest branch represents a buoyant mode, 
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followed by a shear mode at higher va'uea of G. The interface mode is 
shown in more detail in Figure 2, where values of as a function 
of u are also included, *tote that < 0, meaning that the interface 
deformation propagates upward, whereas the base flow is downward near the 
interface. The buoyant mode is shown in Figure 3 and has a a ± value 
three orders of magnitude greater than the interface mode. There is also 
a shear mode with a minimum Grashof number of 6200 shown in Figure 4. 

The variation of the shear and buoyant* modes of instability as a 
function of Prandtl number for the flow between two vertical rigid walls 
has been described by Gershuni and 2hukhovitskii [10]. In Figures 5-7, we 
show the behavior of the various modes as a function of Prandtl number 
for the crystal-melt interface. Figure 5 shows the minimum Grashof number 
(as a function of u) as the Prandtl number varies. In Figures 6 and 7, 
ve have plotted 0 { and w, respectively, as a function of Prandtl 
aumber; Che o± and u curves correspond to minimum values of Che 
Grashof number as a function of u. The buoyanC mode, which has a large 
value of oj_ and a large wave speed (comparable Co the unperturbed 
flow velocity), is hardly affected by the presence of the crystal-melt 
Interface, and is essentially identical to the buoyant mode for a rigid 
boundary. For high values of the Prandtl number the Grashof number for 
the interface mode is two orders of magnitude lower than the Grashof 
number for the shear mode of instability. As the Prandtl number decreases, 
the interface and shear modes approach each other, while the Grashof 
number for the buoyant modes rises abruptly and is much larger than the 
Grashof number for the shear mode for Prandtl numbers less than 10. 
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At a Prandtl number of 2.5, a plot of Grashof numbu. .i. t tiun 

of a), shown In Figure 8, has two minima, which correspond to the two 
Grashof numbers appearing in Figure 3. The size of each of the last 
three terms appearing In eq. (15) are shown in Figure 9; that is, the 
shear term (dashed line), the buoyant term (dot-dashed line), and the 
interface term (solid line). Positive values Indicate a destabilizing 
influence since the term tends to Increase the kinetic energy of the 
disturbance to the main flow, whereas negative values Indicate a 
stabilizing influence. At wavenumbers less than unity, the interface 
term dominates the instability, while for wavenumbers greater than two, 
the shear term is driving the disturbance. 

As the Prandtl number is decreased, Figures 10-17, the second minimum 
in Grashof number, corresponding to the shear mode, disappears (cf. Figure 
3). Both shear and interface' terms are significant for the remaining 
minimum. As P is further reduced, the wavenumber at the minimum Grashof 
number increases and the effect of the Interface term becomes less and 
less significant. For Prandtl number P * .01, Figures 16 and 17, the 
crystal-melt interface term is negligible and the shear mode drives the 
disturbance. 

In Figures 18-20, we show results for lead with L » 0.25 cm, A a 625, 
P - 0.0225, P 8 « 0.012, k r « 0.535, and y " e - 0, for both the 
crystal-melt interface and a rigid boundary. The Grashof number as a 
function of u is essentially independent of the type of boundary. Fo^ a 
rigid boundary ** 0, while for the cryatal-melt interface « -2.7 
at the minimum value of the Grashof number of 7630. This corresponds to a 
period of oscillation of 60 s. 
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Planar Crvstal-Melt Interface 
Succmonitrile (P=22.3) 
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Fig, 1, The Grashof number, G, at the onset of instability as a function 
of the spatial frequency, to, of a sinusoidal peicurbation for 
Prardtl number, P, of succinonit rile , viz., P ■ 22.8. "*he 
three curves correspond to different modes of instability. 
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g. 2. The Grashof number (solJd curve, left ordinate) and the 1 
part of the time constant, o±, (dashed curve, right ordin 
as a function of the spatial frequency, u, of a sinusoids 
perturbation for Prandtl number of 22.8. The solid ".urve 
identical to the lower curve of' Fig. 1. 
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Fig. 4. The Grashof number (solid curve, left ordinate) and the imaginary 
part of the time constant, oj., (dashed curve right ordinate) 
as a function of the spatial frequency, u, of a sinusoidal 
perturbation for Prandtl number of 22.8. The solid curve is 
identical to the top curve of Fig. 1. 
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PLRNRR CRYSTRL-MELT INTERFACE 



Fig. 5* The mi glaum Gras ho f number (as a function of spatial frequency) 

as a function of Prandtl number. For large values of the Prandtl 
number, three different modes of instability occur. 




i 


PLANAR CRYSTAL-MELT INTERFACE 



P 


Fig. 6. The imaginary part of the tine constant, a^, (corresponding 
to the minima values of the Gras ho f number as a function 
of tu) as a function of Pr&ndtl number. 
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POWER 


PLANAR CRYSTAL-MELT INTERFACE 



Fig. 9. The interface (solid curve), shear (dashed curve), and buoyant 
(dot-dashed curve) contributions to the rate of change of the 
average kinetic energy (po^/er) of an instability as a function 
of spatial frequency for Prandtl number of 2.5. 



PLANAR CRYSTAL-MELT INTERFACE 



Fig. 10. The Grashof nuaber, G, at the onset of instability as a function 
of the spatial frequency, w, of a sinusoidal perturbation for 
Prandtl nuaber, P, of 2.0. 
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POWER 


PLANAR CRYSTAL-MELT INTERFACE 



Fig. II. The interface (solid curve), shear (dashed curve), and buoyant 
(dot-dashed curve) contributions to the rate of change of the 
average kinetic energy (power) of an instability as a function 
of spatial frequency for Prandtl number of 2.0. 
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POWER 


PLflNflR CRYSTRL-MELT INTERFACE 



Pig. 13. The interface (solid curve), shear (dashed cvcve), and buoyant 
(dot-dashed curve) contributions to the race of change of the 
average kinetic energy (power) of an instability ? a function 
of spatial frequency for Prandtl number of 1.0. 
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PLANAR CRYSTAL-MELT INTERFACE 



Fig. 16. The Grashof number, G, at the onset of instability as a function 

of the spatial frequency, to, of a sinusoidal perturbation for 

Prandtl number, F, of 0.01. 
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POWER 


PLRNRR CRYSTAL-MELT INTERFACE 



Fig. 17. The interface (solid curve), shear (dashed curve), and buoyant 
(dot-dashed curve) contributions to the rate of change of the 
average kinetic energy (power) of an instability as a function 
of spatial frequency for Prandtl number of 0.01. 
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Fig. 18. The Grashof lumber, G, at the onset of instability as a fund 
of the spatial frequency, «, of a sinusoidal perturbation fo 
Prandtl number, P, of lead, viz., P - 0.0225. The solid cur 
corresponds to a rigid interface vhile the dot-dashed curve 
corresponds to a crystal-melt interface. 
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Planar Oystal— Melt Interface 
Lead (P=0.0225) 



Fig- 19- The Grashof number (solid curve, left ordinate) and the imaginary 
part of the time constant, , (dashed curve, right ordinate) 
as a function of the spatial frequency, w, of a sinusoidal 
perturbation for Prandtl number of 0-0225. The solid curve is 
identical to the dot-nlashed curve of Fig. 18. 
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Task 3 


Measurement of The high Teagjerature Theraophysical 
Properties of Tungsten Group Liquids and Solids 

0. W. Bonnell 

Materials Chemistry Division 
Center for Materials Science 

SUMMARY 


The General Electric Space Systems Division (GE) developed prototype 
"gulp® calorimeter has been tested in joint Rice University (RICE) and The 
National Bureau of Standards (NBS) experiments by dropping levitated tungsten 
solid at 2800 K, and found to function mechanically to perfection. The steps 
necessary to establish this device as a working calorimeter for both ground 
and potential space applications Is discussed. The central problems remaining 
are the addition of calorimetric temperature measurement hardware precise to 
0.001 K and resistant to electromagnetic interference, the development of 
techniques to achieve 0.1 percent enthalpy determination uncertainty with 
expected 18 K temperature rises, and direct calibration of the calorimeter. 

The concentration of effort to measurements at the melting point of 
liquid tungsten have yielded the first reported direct measurement of the heat 
of fusion of tungsten. The value obtained was 53.0 ±2.3 kJ/mole (12.7 ±0.5 
kcal/mole) and, In comparison with estimates and previous indirect methods, is 
considered to be a current best value. 
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INTRODUCTION 


The primary objective of this task is to develop, in conjunction with 
General Electric Space Systems Division (GE) and Rice University (RICE), the 
techniques, methodology, and experimental apparatus for measuring theroophysi - 
cal properties (e.g. heat of fusion, enthalpy increment, heat capacity) of 
-aterials at extremely high temperatures. Electromagnetic levitation of 
conductive samples using auxiliary heating techniques such as electron-beam 
heating has been coupled to a drop-type block calorimeter system. Specific 
application to tungsten has been a prime area of interest because it is the 
highest melting element. Measurements of tungsten are crucial for estimation 
of the theroophysi cal properties of othe elements by correlation. The third 
long period transition elements have all proved to be difficult to measure 
directly, and little data is available for VIB and VI IB elements in spite of 
their importance both scientifically and for engineering purposes. Accurate 
data are necessary to encourage further theoretical investigation of the 
liquid state at high temperature. Modern computer capacities promise new 
models based on detailed evaluation by molecular/atomic dynamics [1], but will 
require validation by experimental determination of definitive therac'hysical 
properties. 

The technical difficulty of this effort has provided a measure of the 
dividing line between experiments which can be performed in the terrestrial 
environment, and those which require the absence of gravity for successful 
accomplishment. Continued advance of ground based efforts seems at present 
the clearest means to successful future experiments in space, as considerable 
hardware development is still necessary to take advantage of the low-gravity 
environment. 
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GULP CALORIMETRY 

In the initial development of the levitation technique, it was recognized 
that the technique was a method to gain the advantage of non-material support 
In a one gravity field, giving many of the advantages of the space micro- 
gravity environment. The technique's application to the measurement of the 
heat content of levitated materials added a separate requirement, whose 
simplest solution is the "drop" calorimeter technique, in which gravity 
provides a rapid and hence approximately isothermal transfer of the sample 
from the high temperature state into contact with a calorimetric mass whose 
temperature is near the reference state. A significant problem in conducting 
similar high temperature calorimetry in space is the transfer of a sample to 
the calorimeter without the aid of gravity and without giving up the advantages 
of the containerless environment. The problem is not insuperable. Acoustic 
positioning techniques are capable of manipulating levitated objects [2-4] but 
the acoustic techniques still appear to have severe materials problems to 
overcome at the very high temperatures which are of greatest interest, and are 
most difficult to achieve in a controlled fashion on the ground. Electrostatic 
techniques are possible [5] and can be used in a vacuum, an important advantage 
where some methods of auxiliary heating need to be used. Less well conceived 
and apparently untested are techniques using an electromagnetic pulse to move 
a sample. Such a method would be difficult to develop in a one-g environment 
in conjunction with RF levitation. Since the design of ground based levitation 
coils is still somewhat of an art, another method was sought. In 1974, this 
investigator made the suggestion that it would be more straightforward to 
attempt to move the calorimeter instead, and J. L. Margrave of Rice University 
coined the name "gulp" calorimeter for this technique. 
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The advent of Space Shuttle experiment opportunities has made this 
problem an important one, as high temperature calorimetry is one of the most 
obvious kinds of science certain to benefit from the advantages of space for 
experimentation. General Electric Space Laboratories has undertaken a ground 
pilot development of a "gulp" calorimeter. NBS expertise has been sought in 
this connection in an advisory capacity, since use of the device on the ground 
would almost certainly occur within the framework of the current GE/RICE/NBS 
levitation effort. The final development by GE [6] is a excellent application 
of mechanical engineering ingenuity, and closely mimics the behavior a space- 
based device must have. Certain compromises were necessary. A primary one is 
the physical size of the calorimeter. Traditional drop calorimeters are 
designed for total heat rises of the order of one degree to minimize the 
temperature difference between calorimeter and surroundings. The mathematical 
basis for data analysis of the thermal behavior of an isoperibol (drop-type) 
calorimeter [7,8] invokes a linear approximation of. the combined effects of 
the exponential Newton's law of cooling (conduction/convection) and the power 
law radiation transfer process. This approximation is generally excellent 
(better than 0.0005 average percent deviation per 10 second measurement period 
over 10-20 minute intervals where calorimeter block/jacket temperature differ- 
ences are 1-2 degrees [8]). Th’is linear relationship is used to determine 
when the calorimeter has equilibrated with the heat input from the sample and, 
more importantly, to correct for the heat contribution to/from the calorimeter 
during the equilibration time of calorimeter with sample. For typical 
equilibration times of ca. 3000 seconds, worst case error from this source is 
less than 0.1 percent and thus negligible in measurements of a few tenths 
percent overall accuracy. 
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The GE calorimeter prototype was designed to be tested on the ground, and 
to be a reasonable model (i.e., physical overall dimensions and system total 
mass) for a flight system. In addition, capture times were required to be as 
short or shorter than ground experimental arrangements, to minimize radiation 
losses from the sample. A pnuematic actuation system was chosen by GE as the 
most reasonable system capable of the required speed of operation with 
significant masses. This choice limited the moving mass of the calorimeter to 
the order of one kg, approximately 5 percent of the current ground 
calorimeter. Predicted heat rises for 10 gm samples of freezing liquid 
tungsten at 1695 K are about 18 K. No isoperibol calorimeter has apparently 
been reported using such large heat increases. By the above analysis, one 
would assume that existing data analysis techniques would cause the error from 
such large calorimeter/jacket temperature differences to contribute at least 
one percent to the uncertainty of the heat measurement. In addition, wide- 
spread practice has determined, for calorimeters with 1 atmosphere gas 
pressure (at 1 gravity) in the jacket region, that a jacket/calorimeter gap of 
the order of 1 cm yields the best uniformity of convection driven temperature 
transfer. This was a difficult structural requirement, and was ignored for 
prototype purposes, as this system is expected to operate in vacuimi. Operated 
in atmosphere the gulp calorimeter system would require careful calibration to 
establish its true heat capacity behavior and the thermal effects of movement 
as a function of differential temperature. Even restriction to operation in 
vacuum below the 10 Pa (0.0001 atm) point, where gas thermal conduction 
becomes much less significant, will still require attention, as the radiation 
transfer, and incidental mechanical connection conduction will not be 
negligible. At the moment, there is just too little experimental data 
available to make such corrections a priori. 
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There are two possible solutions to this problem: 1; new numerical 
techniques for evaluation of the true thermal exchange process; and/or 2) 
controlling the jacket temperature to track the block temperature - a form of 
adiabatic calorimetry. For the first, both the radiation and conduction 
processes have well defined functional forms, and a direct determination of 
the net transfer coefficients for the two processes should be relatively easy 
to model from experimental calibrations, although separating the two processes 
may prove difficult. The hardware solution, using the surface temperature of 
the calorimeter as a control signal to heaters controlling all "jacket" 
surfaces near the calorimeter at tne calorimeter surface temperature is quite 
attractive. This adiabatic wall process is used extensively in calorimetry 
where large temperature excursions are routine (e.g. low temperature calorimetry 
and differential scanning calorimetry) and involves only existing commercially 
available technology. Although the gulp calorimeter system is considerably 
more massive, and has both a much more complicated physical geometry and more 
complicated thermal pathways than the above examples, a combination of the 
above two methods should make high accuracy calorimetry possible with the 
prototype design. 

The GE designed calorimeter uses a pnuematic ram w»th special high speed 
valves to move the calorimeter from its rest position to the capture region. 

For ground based testing, this location is just below the work coil, with 
clearance for the momentum opened calorimeter gate. Timing and capture trials 
with levitation heated samples of tungsten were conducted with the prototype 
instrument as a joint GE/RICF./NBS experimental effort. The calorimeter ram 
and gates were instrumented by GE to allow oscillograms of the time trajectory 
of the calorimeter to be obtained and to detect when, in that trajectory, the 
block radiation gates were fully open [6]. A simple time delay relay was 
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inserted in the signal line from the di'op detector (originally constructed for 
tungsten experiments) which triggers the capture cycle. The delay was adjusted 
to retard the pnuematic impulse about 0.02 sec until the timing oscillogram 
intersected the falling sample location in the middle of the gate open time. 

The block radiation gate is counterweighted to open inertia^'y as. the calori- 
meter accelerates, and to close by momentum as tre caloi ' .r stops at 

maximum excursion. The calorimeter system adapter to tht .tation chamber 
was designed to limit excursion to a location approximately 4 cm below the 
levitation point (about 2 cm below the drop detector location), allowing 
clearance for the coil and pedestal assembly [for experimental details, see 
reference 9]. The calorimeter mount has a pair of overlapping radiation 
shutters which are opened by cams moving with the calorimeter. The calorimeter 
is just beginning to open the block gates as it passes the point of opening 
the main external shutters, and the block gate open time of more than 0.05 sec 
allowed several centimeters leeway in sample location at "capture" (defined 
for this purpose as the time the sample crossed the gate hinge plane and would 
be retained by the closing gates, rather than possibly having some of the 
sample ejected). A capture point corresponding to a drop distance of about 7 
cm (0.12 sec) was selected. Actual testing with a hot levitated sample of 
tungsten (no e-beam heating) was successful. This system halves the drop time 
of the current system, and would presumably improve drop success significantly. 
It is still to be determined how robust the system is with missed drops, and 
how well the timing repeats as the system wears. 

The remaining development needed with this new system is the actual block 
temperature determination. Discussions with GE led originally to the selection 
of a miniaturized commercial 100 chm platinum resistance thermometer. It 
appeared that a custom version of the Hewlett/Packard quartz oscillator based 
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thermometer used for the existing system would have been required to meet 
geometry constraints, allow passage through a vacuum wall, and accommodate the 
calorimeter* movement. The installation used. for the existing RICE/GE system 
is unique to a nonmoving system, and would have been extremely difficult to 
adapt to the gulp calorimeter (requiring a hollow ram and other changes). 

Also, the 10 times larger increase in temperature rise anticipated relaxes the 
measurement constraint by the same factor. Although no commercial sources of 
electronic readout units capable of 0.001 K precision have been ’qcated, it 
should not be necessary to resort to potentiometric techniques to achieve this 
temperature sensitivity. However, a significant problem requiring some 
development exists here, as very small voltagns will be necessary to avoid 
self-heating error and to minimize, input energy to the calorimeter. In 
addition, DC determinations of the resistance of the thermometer are quite 
likely to be seriously disturbed by electromagnetic interference (EMI) for the 
RE. heater. Thus, AC excitation and detection techniques will be necessary to 
determine the resistance of the thermometer. Expected heat inputs for the 
highest' temperature experiments is of the order of 8000 J. Equilibration time 

for the'gulpvcalorimeter will likely be between 1500 and 6000 seconds! .. For 

) 

the upper 1 (mi t estimate, if the thermometer is to contribute less than two 
percent of the total heat to the block, the maximum RMS voltage for excitation 
will be 1.6 volts. This is convenient for modern operational amplifiers, and 
a low, (few 100 Hz) frequency AC nelf-balancing bridge circuit of this sensitiv- 
ity should be within the state-of-the-art. The design may be non-triyial 
because of the electromagnetic interference rejection requirement. Other 
sensors, (e.g. thermocouples or thermistors) have been used, but are generally 
not considered to be as suitable (e.g., because of hysteresis). Temperature 
sensitive integrated circuits (e.g., Analog Devices AD590 or equivalent) might 
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be attractive temperature sensoi-s, given that specialized apparatus to achieve 
the required sensitivity oust be constructed. 

The next steps to complete the development of the “gulp 1 ' calorimeter 
concept can be summarized: 

1) Develop temperature measurement technique to required sensitivity with 
adequate EMI rejection. 

2) Calibrate calorimeter by existing techniques, with specific e^ihasis 
on establishing validity of data analysis and calorimeter performance. 

. 3) Incorporate changes as dictated by the results from 1) and 2), and 
repeat until calorimeter system heat capacity/ temperature response is 
known over operating temperature range to 0.1 % level. 

The unique RICE/GE interaction would seem to be ideal for this task, as 
the capabilities of RICE to accomplish the testing required in 2) appear a 
natural complement to the development expertise of GE in the completion of 
this task. 

»~AT OF FUSION OF TUNGSTEN 

The primary thrust of this project has been to assist RICE and GE in 
developing methodology for high temperature thennophysical property measure- 
ment utilizing the GE coupled RF Levitation/Electron-beam heating facility 
with a RICE supplied and operated calorimeter system. A prime problem in this 
work has been sample temperature determination, as the location of the levita- 
tion process in vacuum, the x-ray hazard posed by electron-beam heating and 
the high vapor pressure of the target material,, tungsten, at its melting point 
all impose stringent requi resents on optical access, and make commercial 
pyrometers impractical. For routine monitoring, GE developed an imaging 
system based on a silicon diode v._dO camera, with selection electronics to 
allow using individual diodes as intensity detectors. Recent efforts [10] 



have established the techniques for converting this type of device to a 

practical pyrometer system, with tracability to the International Practical 

Temperature Scale [11]. The combination of the apparent inability of the 

newly redesigned levitation coil to maintain liquid tungsten, the long cycle 

time of the calorimeter, the inherent noise characteristics of the existing 

camera system, and the difficulty of real time compensation for window obscura- j 

tion have made actual temperature measurements extremely difficult. For that 

reason, and in order to isolate problems in this effort, it was decided to 

concentrate efforts on measurements at the melting point, with the specific 

aim of obtaining the heat of fusion of liquid tungsten. 

Previous reports of this parameter have been primarily based on the 
Tamman Rule [12] (AH /T = 2.3 cal/mole.K {9.6 J/mole.K}, 1 calorie = 4.184 
joule) used in the standard references [e.g., references 13 and 14]. This 
gives the value, 36 kJ/mole, for the heat of fusion using the best estimate 
value for the fusion temperature, 3695 K, determined by Cezairliyan [15]. 

Dikhter and Lebedev [16] reported AH^ of 80 ± 4 cal/gm, corresponding to 62 ± 

3 kJ/mole. This value is probably too high, as it predicts a value for the 
entropy of melting (AS^) of 17 J/mole.K, a value almost 25 percent higher than 
the currently highest known value for a transition metal (molybdenum. 12.9 
J/mole.K. AS o for tungsten is expected to be somewhat higher than for molyb- 
denum, although too much of the data for the transition metals is estimated to 
draw strong conclusions. One would expect, in this region of the transition 
metal portion of the periodic table, that the second-to-third long period 

change in AS would be less than a 25 percent increase. Thus, this value is 
m 

viewed as an upper bound. The report by Wouch, et al., [17] is based on the 
application of a cooling model to experimental data with rather large tempera- 
ture error bounds. The model ignores supercooling, and yielded a heat of 
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fusion of 47.7 kJ/mole, with a corresponding entropy of fusion of 12.9 J/mole.K. 
From their reported fit to the model, a change in AH^ of 8 kJ/mole occurs with 
an increase of only 0.5 times the standard deviation of the best fit and ± 32 
kJ/mole at the 2a point is the variance in the fitted AH n . Although somewhat 
insensitive, this model approach has produced the best estimate of Ah^ prior 
to the direct measurements of this task. 

Results of drops estimated to be 50 percent or more molten out of 25 
drops lor this work are given in Table 1. The estimates of percent melt for 
partially fused samples were based on visual inspection. Sectioning is 
planned to refine these values, but some of these samples may be reused in 
subsequent efforts, as partially melted samples are generally more successful. 
The completely melted sample was weighted 5 times the partially melted data in 
recognition of the uncertainty in the partial melt estimates. Two methods 
were used to derive essentially identical values of the heat of fusion. In 
the first method, the difference between radiation corrected [8] drop heat 
(molar basis) and the enthalpy value for the solid, 119.395 kJ/mole, from 
JANAF [18] was corrected to a 100 percent "melted basis, and a weighted average 
taken. In the second method, a two parameter weighted least square analysis 
of corrected drop heat verses fraction melted was made. The heat of fusion 
was then just the slope of this function, and its intercept is the solid 
enthalpy. That intercept was 119.5 ± 2 kJ/mole, in excellent agreement with 
the JANAF [18] value. Because of this agreement, both methods are essentially 
numerically identical, and the value for the heat of fusion obtained by both 
methods was 53.0 ±2.3 kJ/mole. The statistical uncertainty is of the order 
of twice the anticipated error from temperature uncertainty, excess e-beam 
disturbance of the calorimeter, and normal calorimetric errors. 
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Several factors contributed to the relative success of these efforts: 1) 
nearly flawless operation of the new automatic gate system; 2) improved 
initial stability and low outgassing of the specially cleaned ground samples 
[19], and improved heating and holding characteristics of the slightly smaller 
samples (a set of large samples (1.25 cm diameter, 20 gm) were extremely 
stable levitators, but proved impossible to melt, even on the surface); 3) the 
redesigned coil, which improved levitation stability (only in the case of 
sample N did the coil show any evidence of holding a molten sample, ami that 
sample may have completed melting as it was being lowered in the coil prepara- 
tory to dropping in that case). 

However, a number of problems, related to the age of the GE equipment, 
and the large number of hand operations necessary to ensure successful opera- 
tion resulted in a significant number of failures (at least 1 in 4) due to 
system error, either outright hardware mi sf unctions or human error. The 
remainder of failures to melt are attributable to electron-beam adjustment 
failure. This appears to be the current primary difficulty to obtaining 
successful melting. A number of tests and setup techniques were employed in 
attempts to correlate tuning technique with successful melting. A special 
fixture to simulate high temperature emission grounding was constructed. One 
of the larger balls was mounted oh a 0.32 un diameter tungsten sting, attached 
to a baseplate flange, which allowed the placing of a specimen in the coil for 
beam focusing. This established that the levitation field did not seriously 
disturb beam focus, but pointed up a number of equipment improvements 
necessary to make tuning more reproducible. 

It appears that successful melting occurs only when the electron beam is 
focused so tightly that it can actually “drill" a hole through the sample (see 


figure la and lb). The local heating is apparently considerable, and signi- 
ficant amounts of material are ejected, often with accompanying ar~-induced 
shut-off of the RF generator and loss of sample. It would seem that this 
tight focus and consequent undesirable local heating should not be necessary 
for melting on the basis of estimates of the secondary emission characteristics 
at the melting temperature of tungsten. However .^Tfebility to obs^ve *. 1 
control the beam focus as desired has made systematic testing of focus effects 
difficult. 

A new molybdenum shade ring, to protect the smaller work coil and to 
provide a run-to-run place to tune the electron beam witnout mounting special 
fixtures, will be installed. A number of minor hardware modifications to the 
e-beam control apparatus have also been suggested. These improvements should 
be completed about the time GE finishes the development of a new two-color 
pyroaetry camera system, implementing part of the design proposed earlier 
[19]. In order to improve performance, a gg# series of smaller samples will 
be prepared with a weight of approximately 6-7 gm. In addition, a new coil 
with smaller (1.4 verses 1.6 cm) inside diameter and solenoid geometry to 
maximize the field gradient for these smaller samples and improve liquid 
holding capability has been designed and recommended to GE. 

CONCLUSIONS 

A prototype u gulp u calorimeter system has been developed and tested for 
mechanical performance. The design seems suitable for a potential flight 

experiment series, but a number of developments are still necessary before the 

•««* 

unit operates as a true calorimeter. 

The ground-based effort to determine thermophysical properties for 
tungsten have yielded a value of 53.0 ± 2.3 kJ/mole for the heat of fusion of 
tungsten. This value needs additional experimentation for confirmation, but 
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appears at present to represent a best value for that parameter. Experimental 
difficulties are still a serious problem. It is clear that the electron beam 
performance is the most important obstacle to successful levitation melting. 
Focus, beam shape, and positioning are all quite critical factors, even though 
the more than 3000 watts available is about 7 times the necessary power to 
melt the tungsten samples if absorbed. It is recoumsereHH that GE devote some 
effort to improving this aspect of their experimental facility, with a view 
toward determining the critical factors in electron beam melting of levitated 
specimens. Improvements have been suggested to Improve control and allow some 
experimental tests to be performed. 

It appears that this research still requires some developmental effort to 
achieve its goal of obtaining tungsten liquid enthalpy increments. The 
problems are understood, but will requi re laboratory time to resolve. 
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Table 1. 


Sample 

Weight, gm 

Drop enthalpy, 3 kJ/H 

X melted (est) 

AHm, b kJ/M 

A 

7.895 

158.36 

75 

51.9 

D 

9.740 

148.62 

60 

48.5 

F 

9.473 

148.44 

50 

58.2 

6 

9.596 

158.83 

75 

52.7 

N 

9.468 

172.58 

100 

53.18 

R1A 

9.3249 

159.35 

7b 

53.1 

53.0 + 2. 3 


a Radiation loss corrected 
^computed on 100 X melt basis 
c weighted average/regression - see text 
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la. Run P 



lb. Run T 

Figure 1. Recovered tungsten samples, showing electron beam "drilling" 
of stably levitated samples. Note the obvious signs of 
ejecta surrounding the site of melting. The specimens are 1 cm 
in diameter, and the hole in run T is about 0.7 cm deep. 
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